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Abstract 


The primary goal of kesterite alloying is to allow for fine tweaking of the material's 
characteristics for advanced device engineering. Additionally, it is seen as a viable solution to 
inherent kesterite absorber difficulties such as the Cu/Zn disorder or Sn multivalency. The most 
interesting alloying elements for kesterite are Ag replacing Cu, Cd replacing Zn, and Ge 
replacing Sn for cationic substitution, as well as Se replacing S for anionic substitution. This 
research work investigates the effect of alloying CZTS with Silver (Ag) (Cation) and Selenium 
(Se) (Anion) theoretically using Density Functional Theory (DFT). The compounds were found 
to exhibit indirect bandgap characteristics, with conduction band minima (CBM) and valence 
band maxima (VBM) located between the N and gamma points of the Brillouin zone for pure 
kesterite CuzZnSnSyq and between the N and A points for alloyed AgoZnSnSey respectively. The 
bandgap of around 1.22 eV and 0.78 eV were recorded for the pure and alloyed kesterite 
materials. From the obtained results, a shrink in bandgap was observed due to the presence of 
heavy anion (Se) and cation (Ag) alloying at the same time. It was also found that the 
contribution of different atomic orbitals to the formation of the valence and conduction bands 
is approximately identical for pure and alloyed materials. 
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I. INTRODUCTION 


lobally, research into renewable energy sources, 

including photovoltaic (PV) solar cells, is ongoing to 
fulfill growing energy demand and combat climate change. 
The use of chalcogenide materials such as Cu(In,Ga)(S,Se)2 as 
the absorber layer in solar cells has garnered significant 
interest from the scientific community over the last 40 years, 
culminating in a world record efficiency of 23.35% in January, 
2019 [1]. Recently, a family of similar compounds known as 
kesterite materials enabled the accomplishment of efficiencies 
ranging from 10% to 12.6% [2]. This material family is 
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demonstrated to be a good contender for low-cost, abundant- 
on-earth, high-efficiency solar cells. In 2011, the European 
Commission (EC) announced interest in technologies that do 
not require 'Critical Raw Materials' (CRM). The EC produced 
a peer-reviewed list of CRMs in 2017 that included Ga and In 
[3]. 

Historically, due to the structural similarities between 
kesterite and chalcopyrite compounds, the typical device 
construction for Cu(In,Ga)(S,Se)2 (CIGSSe) was immediately 
extended to CZTSSe by simply substituting a p-type CZTSSe 
thin film for the CIGSSe absorber layer. Typically, CZTSSe 
solar cells are manufactured on soda-lime glass (SLG) 
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substrate covered with a sputtered Mo layer to serve as the rear 
metallic contact. Typical sputter-deposited Mo layer thickness 
is around 500 nm up to | um [4]. After that, the kesterite 
absorber is placed on top of the Mo layer. The quaternary 
systems, Cu-Zn-Sn-S and Cu-Zn-Sn-Se, are derived from CIS 
by substitution of two In atoms for one Zn and one Sn atom. 
This is called isoelectronic substitution and it produces a 
material with many of the same properties as the parent 
compound except that it no longer includes the expensive 
element. This new material still maintains a p-type 
conductivity and a high absorption coefficient (Absorption 
coefficient = a > 10+cm, equivalent to 90% of the incident 
light being absorbed to the surface within 100 nm) [5]. The 
sulfide compound has a bandgap of 1.5 eV, and the selenide 
compound has a bandgap of 0.96 eV [6]. The alloying of 
CZTS and CZTSe to form a pentenary system allows for 
tuning of the material’s bandgap based on the sulfur to 
selenium ratio. Another advantage of the similarity between 
this material and CIS is that CZTSSe can be directly 
substituted into the standard thin film device structure. 
Quaternary CZTS with [2-II-IV-VI4 structure can be 
derived from CIS by replacing In** with Zn?* and Sn**. As 
long as the octet rule is satisfied, it is possible to derive several 
ternary/quaternary semiconductor materials, including 
CulnS2, CuGaS2, CuFeS2, CurZnSnS4, CurFeSnSu, 
CuzCdSnS.4, Cur2ZnGeSy and their selenide counterparts [5]. 
It has been shown theoretically and experimentally that 
tetragonal kesterite is the most stable crystal structure for 
CuzZnSn(S,Se)4 although several metastable structures 
including stannite, Zinc blende, wurtzite, wurtzite-kesterite 
and wurtzite-stannite have also been observed. Kesterite 
CZTS has a direct bandgap of 1.5 eV and a high absorption 
coefficient (>10*) which is close to the ideal for a single- 
junction solar cell. The bandgap can also be tuned by altering 
the S/Se ratio between 1.0 eV for pure selenide and 1.5 eV for 
pure sulfide compounds [7]. 
The scope of this research is to theoretically investigate the 


model of CTZS by replacing the copper (Cu) with Silver (Ag) 
for cation alloying, as well as Selenium (Se) replacing Sulphur 
(S) for anionic alloying. It will be limited to calculating the 
electronic property of both the pure Cu2ZnSnS4 and Co- 
Alloyed AgsZnSnSe, Kesterite within the Density functional 
theory (DFT). 


Il. METHODS 


The calculations were performed on a super-cell structure 
relative to the standard conventional unit cell compounds, 
using first principle calculation by Quantum ESPRESSO 
simulation package code within DFT. Within the Perdew- 
Burke-Ernzerhof generalized gradient approximation (PBE- 
GGA) [8], a pseudo potential was utilized to approximate the 
exchange-correlation potential for both pure Cu2ZnSnS, and 
alloyed Ags»ZnSnSe4 Kesterite. The Maxfessel-Paxton 
smearing method, was used for integrals, while the 
Monkhorst-Pack technique was used to integrate the Brillouin 
zone [9], with 10x10x1 and 12x12x1 k-points grids for the 
pure and alloyed compound respectively. 

For alloying in the Copper atomic site, all the copper atoms 
in the cells were replaced with Ag as cation alloyed and all the 
sulphur atoms were replaced with Selenium as anion alloyed. 
The cell dimensions are constant throughout the calculations; 
however, the atomic locations are completely loosened for all 
calculations using Broyden-Fletcher-Golfarb-Shanno 
(B.F.G.S.) algorithms, until the forces acting on the atoms 
were less than 0.001 eV [10]. 


Ill. RESULTS AND DISCUSSION 


A. Electronic Band Structure 

Electronic band structure calculated along the high 
symmetric point's [—>H—N—I—A in irreducible Brillouin 
for the pure kesterite Cu2ZnSnSyq and alloyed AgoZnSnSeq are 
shown in Fig. 1, with Fermi level at 0 eV. 


Fig. 1 (a) Calculated Band structure of pure kesterite CuzZnSnSyq (b) of Alloyed AgoZnSnSeq4 
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The values of the respective band gaps are clearly indicated 
in Fig. 1 for each of the related band configurations. The 
graph's entire negative energy region corresponds to the 
valance band, whereas the positive energy portion 
corresponds to the conduction band. Both systems exhibited 
indirect bandgap characteristics, with conduction band 
minima (CBM) and valence band maxima (VBM) located 
between the N and gamma points of the Brillouin zone for 
pure kesterite Cu2ZnSnSy, and between the N and A points for 
alloyed AgsZnSnSeq respectively. In pure kesterite 
CuzZnSnS4, a larger bandgap of around 1.22 eV was 
discovered; this value is consistent with [6] and is similar to 
some other pure kesterite systems as reported by [11], while 
the alloyed system exhibits a relatively smaller bandgap of 
about 0.78 eV as compared to pure kesterite. The observed 
shrink in bandgap due to the presence of heavy anion (Se) and 
cation (Ag) alloying at the same time is because VBM of the 
pure CuzZnSnS, is formed of hybridized Cu 3d states and S p 
states, the presence of Ag-3d and Se p states would 


Dos (arbt. units) 


Energy (ev) 


undoubtedly affect the shrinkage. This demonstrates 
unequivocally that the lighter atomic levels in heavier anions 
and cations contribute to VBM at higher states, hence 
lowering the bandgap [12, 13]. This estimated band gap value 
for the alloyed is found to be smaller than the bandgap values 
empirically observed for certain similar alloyed kesterite. 
However, this is consistent with previously published band 
gap estimates derived from density functional theory. 

To gain a better understanding of the electronic properties, 
the total and partial density of states (DOS) for pure kesterite 
Cu2ZnSnS4 and alloyed AgoZnSnSeq were calculated and 
displayed in Fig. 2 within the energy range of -6 eV to 6 eV, 
thereby validating the calculated electronic band structure. As 
observed in the graph, the contribution of different atomic 
orbitals to the formation of the valence and conduction bands 
is approximately identical for pure and alloyed materials. For 
pure and alloyed materials, the Cu/Ag-d orbital and the S/Se- 
p orbital are the primary contributors to electronic states near 
the Fermi energy of the valance band. 


——Zn-s—— Zn-p—— Zn-d 
AAO A. 


2 Energ (ev) 


Fig. 2 Total and partial densities of states (DOS) were calculated for (a) pure kesterite CuzZnSnS, and (b) alloyed 
AgoZnSnSea. 


While electronic states in the conduction band at the Fermi 
energy are dominated by Sn-p and S/Se-p orbitals coupled 
together, all cations in pure kesterite and alloyed Ag2ZnSnSe4 
have a tetragonal symmetry due to the presence of S/Se atoms. 
Thus, as long established by their distinct DOS spectra, 
electronic states are the significance of hybridization of 
individual cation and anion atomic orbitals [12, 13]. 


IV. CONCLUSION 


Using the first-principles methods, an alloyed compound of 
CZTS has been successfully modeled. The compounds were 
found to show an indirect bandgap characteristic with 
conduction band minima (CBM) and valence band maxima 
(VBM) located between the N and gamma points of the 
Brillouin zone for pure kesterite Cu2ZnSnS,4 and between the 
N and A points for alloyed AgoZnSnSeq respectively. The 
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estimated band gap value for the alloyed was found to be 
smaller than the bandgap values empirically observed for 
certain similar alloyed kesterites. This is observed to be 
consistent with previously published band gap estimates 
derived from density functional theory. 
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